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I present the results of a new approach to the intensity and photon density spectrum of the intergalactic 
background light as a function of redshift using observational data obtained in many different wavelength 
bands from local to deep galaxy surveys. This enables an empirical determination of both the EBL and its 
obscrvationally based uncertainties. Using these results one can place 68% confidence upper and lower limits on 
the opacity of the universe to 7-rays, free of the theoretical assumptions that were needed for past calculations. 
I compare our results with measurements of the extragalactic background light, upper limits obtained from 
observations made by the Fermi Gamma-ray Space Telescope, and new observationally based results from 
Fermi and H.E.S.S. using recent analyses of blazar spectra. 



1. INTRODUCTION 

Past work on estimating the spectral and redshift 
characteristics of the intergalactic photon density, 
generically referred to as the EBL, have depended 
on various assumptions as to the evolution of stellar 
populations and dust absorption in galaxies. A de- 
tailed review of the problem has been given by Dwek 
& Krennrich [lj. There have also been attempts to 
probe the EBL from studies of blazar spectra @, [j| 
an approach originally suggested by Stecker, De Jager 
& Salamon In this paper, I present the results 
of a new, fully empirical approach to calculating the 
EBL, and subsequently, the 7-ray opacity of the Uni- 
verse from pair production interactions of 7-rays with 
the EBL. This approach of M.A. Malkan, S.T. Scully 
and myself, hitherto unattainable, is now enabled by 
very recent data from deep galaxy surveys spanning 
the electromagnetic spectrum from millimeter to UV 
wavelengths and therewith using galaxy luminosity 
functions for rcdshifts up to z = 8. This new ap- 
proach, in addition to being an alternative to the ap- 
proaches mentioned above, is totally model indepen- 
dent; it does not make assumptions as to galaxy dust 
or star formation characteristics or as to blazar emis- 
sion models. It is also an approach uniquely capable 
of delineating empirically based uncertainties on the 
determination of the EBL. More details of the work 
presented here have now been published [5|.(See also 
Ref. @.) 



2. GALAXY EMISSIVITIES AND 
LUMINOSITY DENSITIES 

The observationally determined co-moving radia- 
tion energy density u v {z), is derived from the co- 
moving specific emissivity £ u {z), which, in turn is 
derived from the observed galaxy luminosity func- 
tion (LF) at redshift z. The galaxy luminosity func- 



tion, <J>j,(Z/), is defined as the distribution function of 
galaxy luminosities at a specific frequency or wave- 
length. The specific emissivity at frequency v and 
redshift z (also referred to in the literature as the lu- 
minosity density, PlJ) , is the integral over the lumi- 
nosity function 



dL u ^(L u ;z) 



(1) 



In view of the well known difficulties in predicting 
galaxy LFs based on galaxy formation models Q , our 
working philosophy was to depend only on observa- 
tional determinations of galaxy LFs in deriving spe- 
cific emissivities. 

There are many references in the literature where 
the LF is given and fit to Schechter parameters, but 
where pl v is not given. In those cases, we could not 
determine the covariance of the errors in the Schechter 
parameters used to determine the dominant statistical 
errors in their analyses. Thus, we could not ourselves 
accurately determine the error on the emissivity from 
equation ([T]). We therefore chose to use only the pa- 
pers that gave values for pl,X z ) = £v(z) with errors. 
We did not consider cosmic variance, but this uncer- 
tainly should be minimized since we used data from 
many surveys. 

The co- moving radiation energy density u v (z) is 
the time integral of the co-moving specific emissivity 
£u{z), 

dz'£„,(z')^-(z')e- T ^ z ' z \ (2) 

where v' = v(\ + z')/(l + z) and z max is the redshift 
corresponding to initial galaxy formation [8] , and 



dt 



(z) = [ff (l + zVOv + On(l + 2) 3 



dz 

with tt A = 0.72 and tt m = 0.28. 



(3) 
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The opacity factor for frequencies below the Lyman 
limit is dominated by dust extinction. Since we were 
using actual observations of galaxies rather than mod- 
els, dust absorption is implicitly included. The re- 
maining opacity r v refers to the extinction of ionizing 
photons with frequencies above the rest frame Lyman 
limit of VLyL = 3.29 x 10 15 Hz by interstellar and in- 
tergalactic hydrogen and helium. It has been shown 
that this opacity is very high, corresponding to the ex- 
pectation of very small fraction of ionizing radiation 
in intergalactic space compared with radiation below 
the Lyman limit In fact, the Lyman limit cutoff 

is used as a tool; when galaxies disappear when using 
a filter at a given waveband {e.g., "[/-dropouts", "V- 
dropouts") it is an indication of the redshift of the 
Lyman limit. 

We have therefore replaced equation (|2|) with the 
following expression 



very small fractional wavelength intervals, AA(z). We 
have directly tested this by using numerical trial runs. 



4. THE PHOTON DENSITIES WITH 
EMPIRICAL UNCERTAINTIES 



The 68% confidence band upper and lower limits of 
the EBL were determined from the observational data 
on pl u . We made no assumptions about luminosity 
density evolution. We derived a luminosity confidence 
band in each waveband by using a robust rational fit- 
ting function characterized by 

e n ax + b 

PL » = £ » {Z) = ex* + dx + e (5) 



u v {z) = j*™* dz 1 £„,{*') f z {z')n{v{z') ~ v' LyL ), 

(4) 

where H(x) is the Heavyside step function. 



3. EMPIRICAL SPECIFIC EMISSIVITIES 

We have used the results of many galaxy sur- 
veys to compile a set of luminosity densities (LDs), 
Pl„{z) = £ v {z), at all observed redshifts, and at rest- 
frame wavelengths from the far-ultraviolet, FUV = 
150 nm to the / band, / = 800 nm. The LDs were 
obtained with a wide variety of instruments in many 
different deep fields [5j. Figure 1 shows the redshift 
evolution of the luminosity £ u (z) for the various wave- 
bands based on those published in the literature. In 
order to determine the redshift evolution of the LD 
in each of the bands out to a redshift of ~ 8 where 
only UV data are available, we utilized observed color 
relations to transform data from other bands. We 
have chosen to include all data possible at z > 1.5 in 
order to to fill in the observational gaps for various 
wavebands, mostly at higher redshifts. We used the 
redshift-dependent observations of average galaxy col- 
ors where appropriate in our analysis. In the redshift 
ranges where they overlap, the colored (observational) 
data points shown in Figure 1 for the various wave- 
length bands agree quite well ( within the uncertainties ) 
with the black data points that were extrapolated from 
the shorter wavelength bands using our color relations. 
The observationally determined LDs, combined with 
the color relations, extend our coverage of galaxy pho- 
ton production from the FUV to NIR wavelengths in 
the galaxy rest frame. Our final results are not very 
sensitive to errors in our average color relations be- 
cause the interpolations that we made were only over 



where x = log(l + z) and a,6,c,d,and e are free param- 
eters. 

The 68% confidence band was then computed from 
Monte Carlo simulation by finding 10 5 realizations of 
the data and then fitting the function to the form 
given by eq. ([5]) . In order to best represent the toler- 
ated confidence band, particularly at the highest red- 
shifts, we chose to equally weight all FUV points in 
excess of a redshift of 2. Our goal was not to find the 
best fit to the data, but rather to find the limits toler- 
ated by the current observational data. In order to per- 
form the Monte Carlo analysis of the fitting function, 
a likelihood was determined at each redshift given the 
existing data. The shape of this function was taken to 
be Gaussian (or the sum of Gaussians where multiple 
points exist) for symmetric errors quoted in the liter- 
ature. Where symmetric errors are not quoted it is 
impossible to know what the actual shape of the like- 
lihood function is. We have chosen to utilize a skew 
normal distribution to model asymmetric errors. This 
assumption has very little impact on the determina- 
tion of the confidence bands. The resulting bands are 
shown along with the luminosity density data in Fig- 
ure 1. 

With the confidence bands established, we took the 
upper and lower limits of the bands to be our high 
and low EBL constraints respectively. We then inter- 
polated each of these cases separately between the var- 
ious wavebands to find the upper and lower limit rest 
frame LDs. The calculation was extended to the Ly- 
man limit using the derivative derived from our color 
relationship between the near and far UV bands. The 
co-moving radiation energy density was then deter- 
mined from equation (0} . This result was used as in- 
put for the determination of the optical depth of the 
universe to 7-rays . Our resultant z = EBL is shown 
in Figure 2 and compared with the present data, as 
discussed in Ref. @. 
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Figure 1: The observed specific emissivities in our standard astronomical fiducial wavebands. The lower right panel 
shows all of the observational data used. In the other panels, non-band data have been shifted using observed color 
relations in order to fully determine the specific emissivities in each waveband. The symbol designations are FUV: 
black filled circles, NUV: magenta open circles, U : green filled squares, B: blue open squares, V: brown filled triangles, 
R: orange open triangles, /: yellow open diamonds. Grey shading: derived 68% confidence bands. 



5. THE OPTICAL DEPTH FROM 7-7 
INTERACTIONS WiTH UV-IR PHOTONS 

The photon density 

n(e,z) = u(e,z)/e , (6) 



with e = hv, h being Planck's constant, were calcu- 
lated using equation (0). 

The cross section for photon-photon scattering to 
electron-positron pairs can be calculated using quan- 
tum electrodynamics The threshold for this in- 
teraction is determined from the frame invariance of 
the square of the four-momentum vector that reduces 
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Figure 2: Our empirically-based determination of the EBL together with lower limits and data as described in Ref. Q. 



to the square of the threshold energy, s, required to 
produce twice the electron rest mass in the c.m.s., 



s = 2eE y (l - cos9) = 4m 2 



(7) 



This invariance is known to hold to within one part 
in 10 15 [12J, 113] • With the co-moving energy density 
u u {z) evaluated, the optical depth for 7-rays owing 
to electron-positron pair production interactions with 
photons of the stellar radiation background can be 
determined from the expression [4J 



T(E ,Z e ) 



dz ^- I dx^r I dv (1 + zf 



u dz j,, 

u v (z) 



h v 



<7 77 [s = 2Eohi>x(l + z)], 



where Eq is the observed 7-ray energy at redshift zero, 
v is the frequency at redshift z, z e is the redshift of 
the 7-ray source at emission, x = (1 — cos 9), 
6 being the angle between the 7-ray and the soft back- 
ground photon. 

The pair production cross section er 77 is zero for 
center-of-mass energy yfs < 2m e c 2 , m e being the elec- 
tron mass. Above this threshold, the pair production 
cross section is given by 



It follows from equation (J7J that the pair- 
production cross section energy has a threshold at 
A = 4.75 /im • _E 7 (TeV). Since the maximum A that 
we consider here is in the rest frame I band at 800 
nm at redshift z, and we observe E 1 at redshift 0, 
so that its energy at interaction in the rest frame is 
(1 + z)E~, we then get a conservative upper limit on 
E~ of ~ 200(1 + GeV as the maximum 7-ray 
energy affected by the photon range considered here. 
Allowing for a small error, our opacities are good to 
- 250(1 -I- z)- 1 GeV. The 68% opacity ranges for 
z = 0.1, 0.5, 1, 3 and 5 are plotted in Figure 3. 

The widths of the grey uncertainty ranges in the 
LDs shown in Figure 1 increase towards higher red- 
shifts, especially at the longest rest wavelengths. This 
reflects the decreasing amount of long-wavelength 
data and the corresponding increase in uncertainties 
about the galaxies in those regimes. However, these 
uncertainties do not greatly influence the opacity cal- 
culations. Because of the short time interval of the 
emission from galaxies at high redshifts their photons 
do not contribute greatly to the opacity at lower red- 
shifts. Figure 3 shows that the opacities determined 
for redshifts of 3 and 5 overlap within the uncertain- 
ties. 



cr 77 (s) 



16 



2/3(/3 2 -2) + (3-/3 4 )ln 



1-0 



where or is the Thompson scattering cross section 
and = (1 - 4m 2 c 4 /s) 1/2 PI- 



6. Results and Implications 

We have determined the EBL using local and deep 
galaxy survey data, together with observationally pro- 
duced uncertainties, for wavelengths from 150 nm to 
800 nm and redshifts out to z > 5. We have presented 
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Figure 3: The empirically determined opacities for 
redshifts of 0.1, 0.5, 1, 3, 5. The dashed lines are for 
r = 1 and t = 3 [|. 



our results in terms of 68% confidence band upper and 
lower limits. In Figure 2, we compare our 2 = result 
with both published and preliminary measurements 
and limits. As expected, our z = (EBL) 68% lower 
limits as shown in Figure 2 are higher than those ob- 
tained by galaxy counts alone, since the EBL from 
galaxies is not completely resolved. 

Figure 4 shows our 68% confidence band for r = 1 
on an energy-redshift plot |14| compared with the 
Fermi data on the highest energy photons from extra- 
galactic sources at various redshifts as given in Ref. 
[15I | . It can be seen that none of the photons from 
these sources would be expected to be significantly 
annihilated by pair production interactions with the 
EBL. This point is brought out further in Figure 5. 
This figure compares the 68% confidence band of our 
opacity results with the 95% confidence u ppe r limits 
on the opacity derived for specific blazars (l5j |. 

In a recent publication, the Fermi Collaboration has 
probed for the imprint of the intergalactic radiation 
fields (the EBL) in the 7-ray spectra of blazars 0, 
an approach originally suggested in Ref. [3]. Their 
result appears to be consistent with our results near 
the low opacity end of our uncertainty range. The 
H.E.S.S. group Q has also looked for such an effect 
in the spectra of bright blazars at energies above 100 
GeV. It follows from eq. (O that such air Cerenkov 
telescope studies are sensitive only to interactions of 
7-rays with infrared photons. The H.E.S.S. group has 
recently obtained a value for the z = EBL of 15 ± 
2 stat ± i sys nW nT 




Redshift 



Figure 4: An energy-redshift plot of the 7-ray horizon 
showing our uncertainty band results [f| compared with 
the Fermi plot of their highest energy photons from 
FSRQs (red), BL Lacs (black) and GRBs (blue) vs. 
redshift EH. 



This compares to our value of 17.5 ± 4.9 nW m 2 sr 1 
at 0.9 |um. 

My colleagues M.A. Malkan and S.T. Scully and 
I are presently continuing our studies of the photon 
density spectrum as a function of redshift into the in- 
frared range using surveys from Hubble, Spitzer, Her- 
schel and other sources. Such studies, together with 
ongoing complementary 7-ray observations of extra- 
galactic sources with Fermi and future observations 
using the Cerenkov Telescope Array, which will be 
sensitive to energies above 10 GeV [Tg] , one can look 
forward to a obtaining better understanding of both 
the EBL and other potential aspects of 7-ray propa- 
gation in the Universe, such as those explored in Refs. 

E3- El- 
Results Online 

Our results in numerical form are available at the 
following link: 



http : / / csma31 . csm. jmu. edu/physics/ 



scully/ opacities . html 
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Figure 5: Our opacity results for the redshifts of the blazars indicated [f| compared with 95% confidence opacity upper 
limits (red arrows) and 99% confidence limits (blue arrows) as given by the Fermi analysis fl5l ]. 
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